混凝土重力坝爆炸荷载数值分析及抗爆性能研究

利用LS-DYNA非线性有限元程序,基于Eulerian和Lagrangian耦合的方法,研究了RHT本构模型模拟的混凝土板在爆炸荷载下的动力反应,并且将数值结果与现场实验结果进行比较,由此说明了RHT本构模型模拟爆炸荷载下混凝土动力反应的有效性。研究了2 t TNT炸药在距离坝体上游面10 m不同起爆深度的情况下,有泡沫混凝土保护层和无泡沫混凝土保护层的坝体动力响应及其损伤状况。计算结果表明,无泡沫混凝土保护层时,坝体上游面主要损伤区域位置总是随起爆深度的增加向坝体底部移动;当上游表面有泡沫混凝土保护层时,坝体上游表面的损伤明显变小,下游面的损伤较无保护层情况也明显减小。表明泡沫混凝土能够有效减小混凝土大坝在爆炸荷载下的损伤,在提高混凝土大坝的抗爆性能方面起到很好的保护作用。     

混凝土重力坝爆炸荷载数值分析及抗爆性能研究

利用LS-DYNA非线性有限元程序,基于Eulerian和Lagrangian耦合的方法,研究了RHT本构模型模拟的混凝土板在爆炸荷载下的动力反应,并且将数值结果与现场实验结果进行比较,由此说明了RHT本构模型模拟爆炸荷载下混凝土动力反应的有效性。研究了2tTNT炸药在距离坝体上游面10m不同起爆深度的情况下,有泡沫混凝土保护层和无泡沫混凝土保护层的坝体动力响应及其损伤状况。计算结果表明,无泡沫混凝土保护层时,坝体上游面主要损伤区域位置总是随起爆深度的增加向坝体底部移动;当上游表面有泡沫混凝土保护层时,坝体上游表面的损伤明显变小,下游面的损伤较无保护层情况也明显减小。表明泡沫混凝土能够有效减小混凝土大坝在爆炸荷载下的损伤,在提高混凝土大坝的抗爆性能方面起到很好的保护作用。     

各向异性对半圆形河谷散射影响的数值分析

半圆形河谷场地可构造为包含河谷的广义子结构和具有规则边界的开挖场地两部分，基于土-结构相互作用SSI原理，建立子结构控制方程。利用比例边界有限元SBFEM求解开挖场地动刚度，解析求解各向异性介质自由场qP波波动，将两者代入控制方程，可求得广义结构的动力响应。与文献中各向同性半空间中半圆形河谷在P波入射下的位移结果对比，验证了方法的精度和有效性。进一步分析了椭圆各向异性和非椭圆各向异性对半圆形河谷在qP波入射下位移分布的影响。数值算例显示，介质的各向异性改变了半圆形河谷散射位移的空间分布，增大了水平向位移的峰值；同时，介质的各向异性加剧了入射角对散射波场位移分布的影响。     

接触爆炸荷载作用下溢流坝的抗爆性能

以黄登重力坝的溢流坝为研究背景,考虑混凝土的高应变率效应,运用Lagrange-Euler耦合算法建立大坝-库水-空气-炸药全耦合数值模型,研究溢流坝在接触爆炸荷载作用下的抗爆性能。分析满库与空库时溢流坝在爆炸冲击波作用下的动力响应及损伤程度,并进一步研究满库时大坝在不同炸点的水下接触爆炸荷载作用下的动力响应及损伤分布。研究结果表明,满库时水下爆炸比空库时爆炸的动力响应及损伤程度大得多;溢流坝的抗爆薄弱部位主要集中在溢流道顶部及坝体上游折坡处。研究溢流坝的抗爆性能时应重点研究满库时水下爆炸对大坝的破坏特性。     

缝水压力和裂缝面接触条件对重力坝裂缝扩展影响

基于多边形比例边界有限元模拟了重力坝裂缝扩展过程,综合考虑了坝-库水-地基相互作用、裂缝面接触和缝水压力的影响。其中该接触模型能够有效防止裂缝面的相互嵌入;该缝水模型除能考虑裂缝面开合速度、缝隙宽度和水流形态对缝水压力分布的影响,还可基于连续性条件模拟缝内水体流动和空化的现象。以Koyna坝裂缝扩展为例,研究了裂缝面接触条件和缝水压力对计算结果的影响。结果显示,考虑接触的影响,震损较为严重。而与假定为常水压力分布的情况相比,实际的缝水压力分布形式对结果影响较为复杂。在裂缝张开时,减小了裂缝扩展的可能,而在裂缝闭合时正好与之相反。     

热传导问题的比例边界等几何分析

提出比例边界等几何分析SBIGA(Scaled Boundary IsoGeometric Analysis)方法来求解热传导问题。SBIGA兼具比例边界有限元和等几何分析的优势,特别适用于求解包含无限域和奇异物理场的问题。该方法造型十分方便,在径向具有半解析性质,仅需在计算域边界上用NURBS基函数自然离散,为实现CAD/CAE无缝融合提供了新的途径,大大节约前处理和计算耗时。此外,SBIGA无需进一步与CAD系统数据交换就可以保型细分。三个基准算例证明了其在热传导分析中的有效性。与传统比例边界有限元相比,SBIGA模型消除了几何模型误差,并显示出更高的计算精度和收敛速度。     

基于非重叠Mortar方法的比例边界等几何分析

在比例边界等几何分析的建模剖分过程中会出现交界面网格非匹配现象,给计算分析带来困难。为了能够处理此类问题,提出了基于非重叠Mortar方法的比例边界等几何分析。该方法能在将全域分解为若干子域时,针对每个子域分别建模和剖分网格,交界面网格无需逐点匹配;采用交叉点修正的非均匀有理B-样条基函数构造Lagrange乘子空间,子域交界面连续条件可通过Mortar条件满足;并根据交界面连续条件进行自由度凝聚,得到对称正定的系数矩阵,可直接求解。分片试验、U形结构和半无限空间上的柔性基础等数值算例验证了本文方法的有效性,计算精度满足要求。     

海上风机动冰力模型试验相似律研究

相似律在模型试验中发挥着非常重要的作用,虽然目前静冰力模型试验的研究比较全面,但还没有一个完整适用的动冰力相似律作为动冰力模型试验基础.本文针对冰激结构振动中动冰力模型试验相似律,基于弹性相似与柯西(Cauchy)相似的联合相似准则,提出了两种动冰力相似律.建立并改进了能够反映冰与结构水平相互作用的相互作用系数.以NREL5MW单桩海上风机为例,通过是否考虑截面弯曲刚度相似分别给出了两种海上风机在海冰作用下原型与模型的相似关系.利用ANSYS/LS-DYNA有限元方法对原型和模型风机进行数值模拟并开展海冰作用下结构的动力响应分析.根据计算结果,基于不同冰厚的冰-海上风机相互作用分别对提出的两种动冰力相似律进行验证,并对两种相似律进行对比.结果表明,后者考虑了弯曲刚度相似虽然可以优化管壁厚度比尺,为模型加工提供方便,然而模拟精度不如前者.研究成果可为寒冷海域冰力模型试验的研究提供理论参考.     

各向异性对半圆形河谷散射影响的数值分析

半圆形河谷场地可构造为包含河谷的广义子结构和具有规则边界的开挖场地两部分,基于土-结构相互作用SSI原理,建立子结构控制方程。利用比例边界有限元SBFEM求解开挖场地动刚度,解析求解各向异性介质自由场qP波波动,将两者代入控制方程,可求得广义结构的动力响应。与文献中各向同性半空间中半圆形河谷在P波入射下的位移结果对比,验证了方法的精度和有效性。进一步分析了椭圆各向异性和非椭圆各向异性对半圆形河谷在qP波入射下位移分布的影响。数值算例显示,介质的各向异性改变了半圆形河谷散射位移的空间分布,增大了水平向位移的峰值;同时,介质的各向异性加剧了入射角对散射波场位移分布的影响。     

基于等几何分析的比例边界有限元方法

提出了一种具有比例边界有限元的半解析特性和等几何分析的几何特性的新方法。该新方法是在比例边界有限元框架中用NURBS曲线或曲面精确描述域边界几何形状,同时域边界位移场采用描述几何形状的NURBS形函数等参构造。这种新方法具有比例边界有限元固有的径向解析特性和NURBS的高阶连续性的优点。数值算例显示,与传统的比例边界有限元相比,基于等几何分析的比例边界有限元方法提高了域边界单元和域内应力场的连续性,减少了计算自由度。应用此方法可以用较少的计算自由度获得更高连续阶和更高精度的位移、应力和应变场。